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In a deep-infrared (ergodic) regime, QCD coupled to massive pseudoreal and real quarks are 
described by chiral orthogonal and symplectic ensembles of random matrices. Using this correspon- 
dence, general expressions for the QCD partition functions are derived in terms of microscopically 
rescaled mass variables. In limited cases, correlation functions of Dirac eigenvalues and distribu- 
tions of the smallest Dirac eigenvalue are given as ratios of these partition functions. When all 
masses are degenerate, our results reproduce the known expressions for the partition functions of 
zero-dimensional a models. 



PACS number(s): 05.40.-a, 12.38. Aw, 12.38.Lg 

The low energy dynamics of quantum chromodynamics 
is dictated by confinement of colored particles and spon- 
taneous breakdown of the chiral symmetry H] . While the 
vector part of the flavor symmetry in a vectorial theory 
such as QCD is protected by Vafa-Witten theorem M, 
its axial part is considered to be maximally broken by 
the quark condensate. Accordingly the pattern of global 
symmetry breaking of QCD coupled to n(> 2) massless 
quarks falls onto either of the following three classes HQ : 



N c > 3, fund. 
N c = 2, fund. 
N c > 2, adj. 



SU(n) x SU(n) -> SU(n), (la) 
SU(2n) -» Sp(2n), (lb) 
SU(n) -> SO{n), (lc) 



where 'fund.' and 'adj.' stand for the representations 
of the SU (N c ) gauge group to which quark fields belong. 
These three classes are assigned Dyson indices [3 = 2, 1,4, 
respectively, according to the anti-unitary symmetry of 
associated Euclidean Dirac operators ||. In the vicinity 
of the chiral limit where the quark masses m are much 
smaller than the typical hadronic scale Aqcd, the the- 
ory is effectively described by a non-linear a model over 
the coset manifold, without involving gluons and quarks 
that are confined. The effective Lagrangian consists of 
terms made of the Goldstone pion field U(x) and the 
quark mass matrix, consistent with the flavor symmetry 
(l|) ((]]. Furthermore, in the 'ergodic regime' where the 
linear dimension L of the system is much smaller than 
the pion Compton length 1 / m, ~ 1/ y/mAQCD, only the 
zero mode of U(x) dominates, so that the effective 'finite- 
volume' partition function in the presence of the angle 
simplifies into a finite-dimensional integral: MM 



Z {2) {9;M) = [ dU exp (Re tr e i9/n MuA , (2a) 
Z {1) (6; M) = [ dU exp ( i Re tr e w/n MUJU T ) , (2b) 

JsU(2n) V 2 / 



Z< 4 >(0;Af)= / dU c^(Retre te/ ^ n) MUU T ) . (2c) 
J SUM V ' 



'SUM 

Here the rescaled mass matrices M are 



M = diag(/ii, 
M = diag(^i, 



i 4 



J = 1„£ 

oo, rrii 



■ ■ ,Hn) 

■ ■ ,Hn) 

J, J = 



09 = 2,4), 

J 09=1), 

1 ' 
-1 

0, m — 'EL^rrii : fixed, 



(3) 



(4) 



and S stands for the quark condensate in the chiral limit. 
The integrals are extended from cosets to full unitary 
groups (for /3 = 1,4), and their Haar measures dU are 
normalized such that J dU — 1. After Fourier transfor- 
mation, the partition function in a sector with a topolog- 
ical charge v takes the form: 

(M) = I dU detT" U exp ( Re tr MU*) , (5a) 

JU{n) 

Zi l \M) = ( dU det~ u U exp( i Re tiMUJU T ], (5b) 

JU(2n) V 2 / 

Zf> (M) = I dU det' 2u U exp ( Re tr MUU T ) . (5c) 

JU{n) 

For the P = 4 case, v is understood to be substituted by 
uN c B . With a help of Itzykson-Zuber-like formula, an 
explicit form of the integral, which can be considered as 
a 'matrix Bessel function', was obtained for (3 = 2 [@,[§: 

n-l det M j- 1 J l/+i _ 1 ( Mi ) 
4 2) (Mt, , Mn) = 2 ^ II 



A(ari, . . ,,x n ) = 



(6) 



i>3 



1 



and namely for identical /k,'s 



det J„ +i _i(M). (7) 



On the other hand, due to technical difficulties, only the 
case with completely degenerated masses was worked out 
for (3 = 1 and 4 II: 



1 



(2n - 1)!! -™+i<hi<"-i 



(8) 



Zl 4) {/!,... lf i) = (n- l)!!Pf (A) (n: even) (9a) 



A'- 



(n: odd), (9b) 



fc=-c 



n — 1 n — 1 

VmCaOj ^— <«,j<^— • 



Although the above expressions are primarily valid for 
n > 2, they can be extended to n = 1, as a separate 
consideration for this mass-gapped case leads to [|| : 



Z\ft(j i )=I v (j l ), 



(10) 



irrespective of the values of (3. Once the partition func- 
tions are computed in a closed form as the above, physical 
quantities such as the topological susceptibility 



/2) _ E^ 2 ^(M) 



(ii) 



are expressed in terms of them. 

An important observation made by Verbaarschot and 
collaborators HJl^l is that the finite-volume partition 
functions (^|) can as well be derived from models much 
simpler than QCD, Chiral Random Matrix Ensembles 
(xRMEs, for a recent review see Ref. [pd|): 



Zf\{m}) = J dW e-^ v ^ w ^ JJdet 

i— 1 



?7li 

-VF 1 * mi 



(12) 



where the integrals are over complex, real, and quater- 
nion real N x (N + v) matrices W for (3 = 2, 1, 4, respec- 
tively, and it is understood for = 4 that twofold degen- 
erated eigenvalues in the determinant are only counted 
once. Their proofs consist of the 'color-flavor transfor- 
mation' |Q that converts the integration variables into 
nx n matrices and the saddle point method under which 



N — > co, mi — > 0, /i; = irp(0)mi : fixed. (13) 
Here p(0) stands for the spectral density of the random 



matrix T> — 



o w 
-wt 



p(\) = (ti 5(\-W)), 



(14) 



at the origin. The %RME is motivated by the mi- 
croscopic theory (Euclidean QCD) on a lattice, with a 
crude simplification of replacing matrix elements of the 
anti-Hermitian Dirac operator p — (<9 M + iA^j^ by 
random numbers T> generated according to the weight 
e — ptvV(W w) _ u nc l er this correspondence, the micro- 
scopic limit (|l3|) is equivalent to Leutwyler-Smilga limit 
([I|) , since the size N of the matrix W is interpreted as the 
number of cites L 4 of the lattice on which QCD is dis- 
cretized, and the Dirac spectral density at zero virtuality 
p(0) is related to the quark condensate by Banks-Casher 
relation £ = np{0)/L 4 ||. These xRMEs are tech- 
nically more suited for the computation of microscopic 
spectral correlations of Dirac operators p ~ T>, than the 
zero-dimensional a models (^|) where inevitable introduc- 
tion of a probe quark pair (partial quenching) |fl4 15 
brings forth additional complication. In the chiral limit 
/i = 0, the eigenvalue correlation function s [[i6|-f20| as 
well as the smallest eigenvalue distributions |21|-|24| have 
been computed for xRMEs with all three values of (3. On 
the other hand, in the presence of finite and generic /it's, 
these quantities have so far been analytically computed 
only for the (3 = 2 case ||-|| (For (3 = 4, the smallest 
eigenvalue distribution in the presence of /x's was numer- 
ically computed in Ref. fl29|l , and the partition function 
with four degenerate /x's was analytically treated in Ref. 
Ippfl ). We shall treat the remaining cases, chiral orthogo- 
nal ((3=1) and symplectic ((3 = 4) ensembles with finite 
mass parameters. This Letter is devoted to the com- 
putation of partition functions of these ensembles. As 
corollaries, we provide closed expressions for the small- 
est eigenvalue distribution for (3=1 and odd v, and the 
microscopic eigenvalue correlation functions for (3 = 4. 
The microscopic eigenvalue correlation functions for the 
(3 = 1 case will be presented in a separate publication 
II- 

We start by expressing the partition function (|lj) of 
the xRME in terms of eigenvalues Xi = Af of the positive 
definite matrix W^W (up to a constant independent of 
m and v): 



Z^\{m})= m< WK{m}), (15) 
\=i ' 

■Jo J=1 

n 

w(x; {m}) = e -0V( x ) x f3(»+i)/2-i ^ x + m 2). (16) 



2 



Since the partition function is even under v — > — v, we 
have set v non-negative integer, without loss of general- 
ity. We note that all spectral correlation functions of the 
orthogonal and symplectic ensembles can be constructed 
from the scalar kernel of the unitary ensemble sharing 
the same weight function w(x) [Bgj33|| . Since the scalar 
kernel in the microscopic limit (|13|) is known to be in- 
sensitive to the details of the potential V(x) either in 
the absence p4] , |3l| or in the presence of finite fi's p5| , 
it suffices to concentrate onto Laguerre (chiral Gaussian) 
ensembles, V(x) — x. This choice leads to 



p(X) = -y/2N- A 2 



(17) 



We sketch the skew-orthogonal polynomial method |3(| 
as employed in Ref. 0], to which we leave details. As 
the final result in the microscopic limit is insensitive to 
the parity of N, we consider only even N henceforth. 



A. orthogonal ensemble 



We use the identity 

JV 



A(xi, . . .,X N ) Yl Y\_( x j - XN+i) 
.7=11=1 

a / \ det Ri_i(Xi) 

A{Xi, . . . ,X N+n ) _ l<i,j<N+n J/ 
A(XN+1, ■ ■ ■ ,XN+n) A(xjV+l, ■ • ■ ,XN+n)' 



(18) 



where Ri (x) is an arbitrary monic polynomial of the «-th 
order, and xn+i = —m? < 0. We take {Ri(x)} to be 
skew-orthogonal 

(i?2i, R2j+i)n — — {R2j+i, R2i)r = hiSij , others = 0, 

(19) 

with respect to the antisymmetric product 



(f,9)n= / dxx i e x g(x) / dyy » e «f(y) 
Jo Jo 

-(/-<?)■ (20) 

When eq.(|l8|) is integrated over x\,...,Xn with the 

weight Y\^ = i(e~ Xi x i 2 ) in a cell < Xi < x 2 < ■ ■ ■ < xjf, 
it can be neatly expressed as a Pfaffian, due to the skew 
orthogonality @ @: 

S £ ) ({™}) = ±- 1 r 7 -* J —^— (n:even) (21a) 



A(m 2 ,...,m2) 
F 

~R T 



(nlff'-'^pff It" 



k=0 



R 1 = i?Ar +n -i(-m 2 ), 



hk 

1 < i,j < n. 



Explicit forms of the monic skew-orthogonal polynomials 
and their norms are known as Bffl: 



(2fe)! d r („_i), 9 > 

-&+^4r 1} (2.) 



(22) 



fc fe = 2- 4fc - ,y (2fc)!(2/c + i/)!. 

In the microscopic limit ( |l3] ) with /ij = 2\ / 2Nrrii fixed, 
the sum over the indices k becomes an integral, and La- 
guerre polynomials approach modified Bessel functions: 



L[ a \x) 



I a (2V^kx') (x = 0(-)<o). (23) 



Then the partition function is expressed as 



4 1} (M) 
#>({/*}) 



i=l ' 

Pf(/) 



Pf 



/ r 
-r T 



* a(m?,-.,m^) 



(rt : even) (24a) 
(n : odd), (24b) 



where 



n— 2 



= (-l) 11 ^— 2" _1 (n- 1)! JJ(2ife+l)! (n : even) 



fc=0 



n-2 



-l)-4- l 2 i — ^ l (n - 1)!! J] (2k + 1) 

fe=0 



(n : odd), 







1 < i, j i < n. 



The constant c n is conveniently determined as the above 
by requiring the small-/i behavior be in accord with that 
of the zero-dimensional a model, 



A(m 2 , . . . ,m 2 ) 



(n:odd), (21b) J] 



— (2fc)! 



fc=0 



(2fc + z/)! 



(/i«l). (25) 



where 
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B. symplectic ensemble 



where 



We concentrate on the case with an even n(= 2a) num- 
ber of flavors and pairwise degenerated mass parameters, 
corresponding to adjoint Dirac fcrmions in the QCD con- 
text. 

We use the identity 

N a 

A(xi, . . . , xn) a Y[ Y\_( X J ~ x N+if 
3=1 i=l 



AT+[f]-l 



1 


X\ 


xf ■ 






x 2N+a-l 





1 


2xi ■ 


■ (2N- 


-a- I) 


x 2N+a-2 


1 


x N 


x% ■ 






2N+a-l 
N 





1 


2xn ■ 


• (2N- 


-a -l] 


2N+a-2 
JV 


1 


xn+i 


2 






2N+a-l 
X N+1 


1 


XN+a 


T 2 

■ L N+a 






2N+a-l 
■ L N+a 



Uxi) 



1o{xn) 



> XN+a) 

Ql{Xi) ■■■ Q2N+a-l(xi) 

Q'liXl) ■■■ Q 2 N+a-l( X l) 



Qi(x N ) 
Q[(x N ) 



]o(xn+i) Qi{xn+i) 
)o(x N+a ) Qi(x N+a ) 



Q2N+h-i(xn) 
Q2N+a-l( x N) 
Q 2 N+a~l(xN+l) 



hN+a-l(xN+a) 



A(x N+1 , . . . , X N+a ) 



, (26) 



where Qi (x) is an arbitrary monic polynomial of the i-th 
order, and x^+i = — mf < 0. We take {Qi(x)} to be 
skew-orthogonal 

{Q%i,Q2j+i)Q = -{Q2j+i,Q%i)Q = Kkj, others = 0, 

(27) 

with respect to the antisymmetric product 



(/,<?>< 



dxx 2 » +1 e~ ix (f(x)g'(x) - f(x)g(x)). (28) 



When eq.(|26D is integrated over x\,...,xn with the 
weight Yl i= i(e~ 4xi xf u+1 ) , it can be neatly expressed as 
a Pfaffian, due to the skew orthogonality (|27|) e3|: 



llS^M Pf(F) 



A(m 2 , . 



(nS IfM 



hA Pf 



F Q 
-Q T 



A(mf,...,m2) 



(a : even) (29a) 
(a : odd), (29b) 



E 

fc=0 



Q2k(-mf)Q 2 k+i{-m 2 J ) - (i <->j) 



= <32N+a-i (-to?), l<«,j<a- 

Explicit forms of the monic skew-orthogonal polynomials 
and their norms are known as |37||: 



>.2h 



(x) 



k\(k + u)\ ^ (2Z-1)!! A2v) 



2 2k 



2^{l + v)\ L21 {AX) 



2 4fe+2 ^2/c+ll 4 ^ 



Q2k+i{x) = - \, k J L^iAx), CUD 
ft* = 2- 8fc - 4 ^- 4 (2fc + l)!(2fc + 2v + 1)!. 



In the microscopic limit with /i^ = 2y2N rrii fixed, the 
sums over the indices k and I become integrals, and La- 
guerre polynomials approach modified Bessel functions. 
Then the partition function is expressed as 



pf(/) 



A( M ?,---,^) 



pf 



/ q 
-q T 

a A(^,...,^) 



(a : even) (31a) 
(a : odd), (31b) 



where 



c a = (-l)^ 1 l[(2k + l) 



k=0 



f 



o M 



,2y 



,2t/ 



2// 



1 < i, 3 < a. 



The constant c a is conveniently determined as the above 
by requiring the small-/i behavior be in accord with that 
of the zero-dimensional a model, 



4 4 v,..., M )^ 2 -n 



2a 



k=0 



(2fc+l)! 
(2k + 2v + l) 



( M «l). (32) 



It remains to confirm whether the above expressions 
for the xRME partition functions agree with those of the 
zero-dimensional a models. By taking all /x's to be iden- 
tical, we obtain for j3 = 1: 



4 



Pf(/) 



Pf 



(n/2-l)n 



-f T 



(n-l) 2 /2 



(n : even) (33a) 



odd), (33b) 



where 



Cn = (- 



fc=0 



n(^-l) („-l)(3r, + ll) rifc-0 + "-)' / J J \ 

(-1) 2 2 s 1 x *~^ (n : odd), 



IUo(2fc + l)' 



/« - / dtt i+j+2 I v+ ^ 1 {t l i)I v+j {t l i) - j), 
Jo 



and for /3 = 4: 



0<*,i<n-l, 



Pf(/) 



2d 



Pf 



/ 9 



-g T 



(a-l)a/2 



where 



/' 



( _ 1)fl ^g(2 fc + l)l 



even) (34a) 



(a : odd), (34b) 



fe=0 



fc! 



9 = 



l ,i 

dtf +J+1 h v+i (2tfi) / duu 3 hu+j{2tu^) 
o Jo 

i 

dtt i I 2 u+ l {2t^), 0<i,j<a-l. 



We have numerically checked that, despite the appear- 
ances, the above expressions are identical to eqs.(^|) and 
(pa|). Together with the [3 — 2 case that has previously 
been confirmed they explicitly show the equiva- 

lence between the xRMEs and the a models in Leutwyler- 
Smilga limit. 



The integral domain in the numerator can be traded 
to [0, oo) with the weight function shifted by s, w(x + 
s;{m}). If the exponent in the weight function + 
1) — 1 is an integer (excluding the case with [3 = 1 and v 
even), we can utilize the 'flavor-topology duality' p9[ 



1 



,m n , 0, . . . ,0 ), (36) 



(H-l)-l 



to express E$f\s;{m}) in terms of the partition func- 
tions: 



f(*+i)-i 



s |li(v /s + m i' • • • > \/ s + m li Vs,...,Vs) 
(roi,...,m n J 



(37) 



Now we change the picture back from Laguerre to chi- 
ral Gaussian, and take the microscopic limit with £ = 
TTp(0)y/s = 2y/2Ns and fa = 2^/2Nrrii fixed. The distri- 
bution of the smallest eigenvalue of chiral random matri- 
ces is then given by the first ^-derivative of E„ 



(0). 



d_ 



e -(/3/8)C 2 x 



if^iVe+A, . . . , 7^+74, . c, . . . , c ) i 
^ ; ■ ( 38 ^ 

& '[pi, . . . ,Hn) ) 

For (3=1 and v odd, eqs.(|3^) and ( |24| ) suffice to express 
the smallest eigenvalue distribution in a closed form. 
This prediction should in future be put in comparison 
with lattice QCD simulations with overlap dynamical 
quarks. 

On the other hand, for (3 = 4, the partition function 
in the numerator falls out of the range of this Letter, 
as the number of additional flavors is odd. A different 
formalism based on Fredholm determinant |l0| might be 
needed in order to overcome this limitation. 



D. correlation function 



C. smallest eigenvalue distribution 

The probability of finding no eigenvalue in the interval 
< x < s is given by 

" frnU(dx jW (xf,{rn}))\A({xW 

(35) 



In the case of even /3, a p-level correlation function 

p 

p(X u . . . , A p ; M) = (H tr5(X k - iV)) (39) 
fc=i 

is expressed by construction as a ratio of partition func- 
tions with n and n + (3p flavors |3^,^l|. After taking 
the microscopic limit, In the = 4 case with a pairs of 
doubly degenerated masses, it reads: 



5 



2n2 \ Z^(jj,i,fii, . . . ,/i a , /i a , . . . , ifi, . . . , . . . , z£ p ) 



fc=l \ i=l 



(40) 



As our derivation of the partition function ( |3l| ) is valid as [17] 



well for negative values of or nf, the above relation- 
ship suffices to express any p-level correlation function in 
a closed form, up to an overall constant independent of 
/i's. 
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